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Abstract 

Attacks by wolves (Canis lupus) on dogs (C. familiaris) presumably are motivated both by preying and elimination of 
potential competitors. Regardless of these alternative motivations in wolves, the risk of attacks might be higher when the 
density of primary prey is low. We examined how many dogs do territorial wolves in Finland kill in relation to the popula- 
tion density of the most abundant ungulates, moose (Alces alces), white-tailed deer (Odocoileus virginianus), and roe deer 
(Capreolus capreolus). Most attacks by wolves on dogs take place in hunting with dogs. The number of wolf-killed dogs 
was in highly significant negative relationship to the population density of white-tailed deer and to total ungulate biomass 
per unit area which is largely determined by the density of white-tailed deer. Our results indicate that abundant wild prey 
would decrease the risk at which wolves attack dogs. On the other side of the coin prevail two hard facts which wildlife 
managers had to take a notice. White-tailed deer, although a potential mitigator of wolf—human conflict, is an alien species 
and a partner in > 6000 traffic collisions annually in Finland. One factor that seemed to increase the risk of wolf attacks on 
dogs is the low ungulate density in regions where moose is the only remarkable ungulate prey. Higher moose densities could 
decrease the risk of attacks, but on the other hand, higher densities could increase the risk of serious traffic collisions and 
browsing damages in forests. 


Keywords Wolf—human conflict - Wolf - Domestic dogs - Ungulate density - Moose - White-tailed deer - Roe deer 


Introduction 2021). Wolf (Canis lupus) attacks on domestic dogs (C. 
familiaris) is a small-scale human-large carnivore conflict 
Attacks on domestic animals have stayed as one of the big- as compared to a conflict owing depredation on sheep and 


gest global challenges to the successful conservation of large other livestock, but a remarkable conflict in some regions 
carnivore populations (Linnell et al. 2001; Cervasi et al. (Fritts and Paul 1989; Kojola and Kuittinen 2002; Backeryd 
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2007; Olson et al. 2015; Bassi et al. 2021; Tikkunen and 
Kojola 2020; Trbojevic et al. 2020). 

The biological environment may affect the risk of 
carnivore-caused damages due to the potential relation- 
ship between the abundance of wild prey population and 
the magnitude of damages wolves might cause by killing 
domestic animals. Prey scarcity hypothesis (Nelson et al. 
2016) assumes that livestock depredations by wolves may 
increase when wild prey exist at low densities in wolves’ 
home ranges. Some studies indicate that high rates of 
livestock depredations are associated with low density of 
wild ungulates (Meriggi et al. 1996; Cervasi et al. 2014; 
Khorozyan et al. 2015; Janeiro-Oteroa et al. 2020), thus 
providing empirical evidence that appears to support prey 
scarcity hypothesis. 

Domestic dogs descend from the wolf and may breed with 
wolves (e.g., Salvatori et al. 2020; Harmoinen et al. 2021). 
Far more often, dog-wolf encounters end up wolves’ fatal 
attacks on dogs (Fritts and Paul 1989; Kojola & Kuittinen 
2002; Kojola et al. 2004; Backeryd 2007; Bassi et al. 2021; 
Tikkunen and Kojola 2020). Domestic dogs differ from live- 
stock animals both in ecological and social dimensions of 
wolf—human conflict. Dogs share the same carnivore guild 
with wolves, and therefore, their principal motivation to attack 
on dogs is not predation for food but the elimination of poten- 
tial competitors although wolves usually also eat from the dog 
they have killed (Fritts and Paul 1989; Kojola and Kuittinen 
2002; Backeryd 2007; Bassi et al. 2021). Occasionally, dogs 
might even constitute a remarkable portion of wolves’ diet 
(Trbojevic et al. 2020). Domestic dogs, irrespective of their 
concrete service roles (hunting, rescue operations, livestock 
guarding, etc.), are usually positioned as family members and 
therefore often regarded as priceless by their owners (Bisi 
et al. 2010). 

Very little is known about potential influences of prey 
abundance on the risk of wolf attacks on domestic dogs. 


Fig. 1 Wolf-killed dogs (black 
dots) and territories of wolves 
in Finland outside the reindeer 
husbandry region in 2016-2020 
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Results from Estonia and eastern Finland provided evi- 
dence that risk is highest in years when prey densities are 
low (Kojola et al. 2022). Herein, we evaluate the prey scar- 
city hypothesis at the wolf territory level in Finland. The 
examination of incidences at the territory level potentially 
provides us with a lot more detailed picture about the effects 
of prey abundancies, and in the current study, we can also 
take into account the estimated number of wolves occupying 
the territory. 


Study area 


Our study area are wolf territories located in Finland out- 
side the reindeer husbandry region from winter 2016/2017 
to winter 2019/2020 (Fig. 1). Predominant landscape is a 
coniferous boreal forest where the primary tree species are 
Scots pine (Pinus sylvestris) and Norway spruce (Picea 
abies). Topography outside the reindeer husbandry region 
is relatively flat; elevation ranges between 0 and 354 m 
above sea level. Of wolves’ most numerous prey, moose 
(Alces alces) exist throughout the study area while white- 
tailed deer (Odocoileus virginianus) only in the southwest. 
Roe deer (Capreolus capreolus) distribution is wider than 
that of white-tailed deer, but population densities are low 
in the east. 


Data 
Wolves 


Wolf territory boundaries (Fig. 1) were assessed from point 
observations for packs and twosomes (pairs) recorded by 
about 2000 voluntary large carnivore contact persons and 
locations by GPS-collared territorial wolves, and non-invasive 
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genetics were used for separation of neighboring packs and 
pairs from each other (Kojola et al. 2018; Heikkinen et al. 
2020). The estimated wolf population size in March increased 
from 150—165 (2017) to 216-246 (2020) during our study 
period (Heikkinen et. 2018, Kojola et al. 2018; Heikkinen 
et al. 2019, 2020). 

Genetic monitoring was performed annually between 56 
and 76% of all territories in 2017-2020 to distinguish neigh- 
boring territories from each other (Heikkinen et al. 2018, 
2019, 2020; Kojola et al. 2018). 

Territory delineations were multiple convex polygons 
(MCP) because only GPS locations would provide data con- 
venient for more advanced home range estimates. The esti- 
mates for pack size are November estimates based on both 
genetic analyses from scats collected by Luke’s field tech- 
nicians and volunteers and the maximal number of wolves 
in point observations made during winter before the end 
of February. The known mortality was taken into account 
in pack size estimation. Genetic analyses and observational 
data provide highly correlated estimates for the pack size 
(Kojola et al. 2018). 


Dogs killed by wolves 


We considered casualties taken place between June 1 and 
March 31 in 2016-2020. We ignored a few cases that took 
place in April and May because most born in the previous 
year depart their natal territory by the end of April (Kojola 
et al. 2006), and thereby, the number of wolves occupying 
the territory during these months was not possible to esti- 
mate. The conclusion is based on movement data for 29 pups 
collared in winter (Kojola et al. 2006). Female wolf gives 
birth to new pups in late April and early May (Kaartinen 
et al. 2010) which results in the formation of a family pack. 
Only the damages accepted in a formal inspection by local 
officials were included in our study. The cases were stored 
in the official register (“Riistavahinkorekisteri’; the register 
of wildlife damages) with dates and geographic coordinates 
(Fig. 1). When we related damages to prey abundance within 
wolf territories, we accepted only the cases that took place 
within the delineated wolf territories (1 =91). These cases 
comprised of 58% of all casualties (n = 156). Nearly all fatal 
attacks on dogs by wolves (97%) took place during autumn 
hunting seasons (2018). 


Abundance of wild prey 


Population densities of moose, white-tailed deer, and roe 
deer were autumn post-hunt estimates by hunting clubs. 
Prey biomass available for the wolves per unit area were 
estimated based on fawn/adult (deer) ratios in hunting bag 
and mean dressed weights of adults and calves or fawns in 
autumn hunting seasons. The resultant biomass estimates 
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were coarse because the age structure in the living popula- 
tion is not known. Because wolf territories were manually 
delineated MCPs, we used biomass estimates only per areal 
unit while not for the total area of the territory. 

The population density of moose varied less among the 
wolf territories than did the density of roe deer especially 
the density of white-tailed deer. White-tailed deer popula- 
tion was highly concentrated in the southwestern wolf ter- 
ritories (Fig. 2). Moose existed in almost all study territories 
(Fig. 2). White-tailed deer was highly concentrated in south- 
western Finland while roe deer was more evenly distributed 
in western Finland (Fig. 2). Only moose population density 
was normally distributed among wolf territories (Fig. 2). 
The mean density of moose was ca. 3.0 moose per 1000 ha. 
The density distribution of white-tailed deer and roe deer 
were biased (Fig. 2) because these deer were absent in most 
territories in eastern Finland (Fig. 2). 


Dogs in wolf territories 


The casualties were tightly connected to hunting seasons. 
The number of dogs within wolf territories was not known. 
The registry for wildlife damages does not provide the 
breed of wolf-killed dog or the personal information about 
the owner of the dog, but an inquiry to communal authori- 
ties provided us with details for 32 cases taken place in 
2017-2020. Only one of these dogs was killed by a wolf 
in house-yard and others in hunting with the dogs. In 20 
casualties (63.2%), the wolf-killed dog was a breed used in 
moose hunting. Local hunting clubs report the number of 
moose-hunting dogs used in moose hunting per season (from 
Oct. 9 through Jan. 15). 


Statistical analysis 


A generalized linear mixed model with a negative binomial 
distribution assumption and log-link function was con- 
structed to model the count response (number of wolf-killed 
dogs in the territory i). The hierarchy of the model consisted 
of two levels: wolves’ territory as the highest level and years 
nested within the territories. The number of years varied 
by the territories, and no serious temporal autocorrelations 
were found between the adjacent years. Thus, the residuals 
could be treated as uncorrelated errors instead of AR(1). The 
model was computed using R-package glImmTMB (Brooks 
et al. 2017). 

In the R function glmmTMB, two variance func- 
tions were tried: NB1 variance = uw(1+a) and NB2 
variance = w(1 + fl (Brooks et al. 2017). NB1-parametri- 
zation suggests the linear mean-variance relationship, and 
NB2-parametrization suggests the quadratic relationship. 
Both parametrizations gave almost the same AIC values 
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White-tailed deer 


European Journal of Wildlife Research (2023) 69:26 


Moose White-tailed deer Roe deer 

20- 20K 20+ 
ra 
3S 
je) 
Oo 

10- 10K 10+ 

0 | L L | ! 0 L ! 0 I I 
0 1 2 3 4 5 6 7 0 10 20 30 40 50 60 70 0 5 10 15 20 


Population density (ind./1 000 ha) 


Fig. 2 Population densities of primary prey species (ind./1000 ha) and wolf territory boundaries A and the distribution of population densities of 
prey species within wolf territories in 2020 B, Finland outside the reindeer husbandry region 


in the models. The NB2-parametrization was used in the 
models, and the dispersion parameters of NB2 models were 
used in the simulations of model fit, giving directly the 0 
-parameter for the simulations using R-function dnbinom. 

R-package glmmTMB allowed also the zero-inflation 
modeling, but the zero-inflation was not an obvious problem 
in the data. It was tested using R-package DHRMa (Hartig 
2021). The test indicated that negative binomial distribu- 
tion worked without any zero-inflation coefficient in all 
the models. The p-values of zero-inflation tests (function 
testZeroInflation) remained far higher from the limit of sig- 
nificance at a 5% risk level. The test compares the number 
of simulated zeroes (function simulateResiduals) against the 
observed ones. 

Variance inflation factor (VIF) was used to check a possible 
multicollinearity of the main effects of the explanatory vari- 
ables. The VIF was computed using R package “performance” 
(Liidecke et al. 2021) and its function check_collinearity. The 
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predictions with standard errors were computed using R-package 
ggeffects (Liidecke 2021) and the effects plots using R-package 
ggplot2 (Wickham 2016). 

All candidate models included the estimated number of 
wolves occupying the territory and the number of moose- 
hunting dogs per area of the territory. Moose-hunting dogs 
are the most common breeds among attacked dogs during 
hunting (Kojola and Kuittinen 2002). The model-specific 
variable for prey abundance was (1) the population density 
of moose, (2) population density of white-tailed deer, (3) 
population density of roe deer, (4) common density of these 
three ungulates and (5) ungulate biomass per unit area, and 
(6) population densities of all these 3 ungulates. All can- 
didate models addressed the status of the territory (pair/ 
pack), pack size estimate (2 for pairs), and the number of 
moose-hunting dogs per unit area (“dog index’). Because 
the area of the territory was not a significant variable, it was 
excluded from the final models We evaluated model fit by 
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using the Akaike information criterion (AIC) (Burnham and 
Anderson 2002). Models with a AAIC <2 were considered 
equally competitive. 


Results 


Of six candidate models (Table 1), the annual number of 
wolf-killed dogs was significantly related to the number of 
moose-hunting dogs per area of the territory (“dog index”) 


Table 1 Parameter estimates 
and tests for the generalized 
linear mixed models using 
negative binomial distribution 


Variable 


assumption for the number of Intercept 

killed dogs (killed by wolves), Wolves, Ind 

when dog index was included in Moose density, Ind./1000 ha 
the models as a covariate. VIF Dog index 


is the variance inflation factor : . 
Variance of territory 


Version 2 (AIC = 327.3, R°=38.4) 


Intercept 
Wolves, Ind 


White-tailed deer density, Ind./1000 ha 


Dog index 


Variance of territory 


Version 3 (AIC = 339.6, R*= 19.3%) 


Intercept 

Wolves, Ind 

Roe deer density, Ind./1000 ha 
Dog index 


Variance of territory 


Version 4 (AIC = 330.5, R?= 30.2) 


Intercept 
Wolves, Ind 


Total density of ungulates, Ind./1000 ha 


Dog index 


Variance of territory 


Version 5 (AIC =331.7, R?=28.1) 


Intercept 
Wolves, Ind 


Total biomass of ungulates, kg/1000 ha 


Dog index 


Variance of territory 


Version 6. (AIC=355.3, R?=47.7) 


Intercept 
Wolves, Ind 
Moose density, Ind./1000 ha 


White-tailed deer density, Ind./1000 ha 


Roe deer density, Ind./1000 ha 
Dog index 


Variance of territory 


Version 1 (AIC =345.0, R?= 12.7%) 
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in two models: the model where the pack size was treated 
as independent variable together with the white-tailed deer 
density or with densities of all three ungulates (Table 1, 
Fig. 3). The number of wolf-killed dogs was related posi- 
tively to the pack size in all six models. 

The model where white-tailed deer density was the 
only variable for prey abundance fitted best with data 
(Fig. 3, model 2 in Table 1). The model where the density 
of white-tailed deer was recovered by the total ungulate 
density (model 4 in Table 1) had only a slightly poorer 


Estimate Std. error z-value p-value VIF 

—3.320E—01 5.369E—01 —0.618 0.536 - 
1.137E-O1 5.599E — 02 2.031 0.042 1.03 

—2.142E—-01 1.297E-01 —1.652 0.098 1.03 

—6.980E—0O1 4.993E—01 —1.398 0.162 1.04 
0.523 

—4.952E—01 3.420E—O1 —-1448 = 0.148 - 
1.312E—01 5.284E— 02 2.483 0.013 1.03 

—8.635E—02 2.291E—02 -3.769 0.000 1.02 

—9.896E-Ol1 4481E—01 —2.208 0.027 1.03 
0.299 

—5.565E—01 3.740E—01 -1488 0.137 - 
1.149E—-01 5.422E—02 2.119 0.034 1.01 

-—1.277E-O1 4.710E—-02 —-2.710 = 0.007 1.00 

—7.246E-0O1 4.608E—01 -1.572 0.116 1.01 
0.420 

—3.444E—-01 3.568E—01 —0.965 0.334 - 
1.224E—01 5.307E—02 2.306 0.021 1.01 

—5.241E—02 1.437E-—02 -3.646 0.000 1.01 

-8499E-O1 4511E—01 —1.884 0.060 1.01 
0.335 

—1.378E-01 3.842E—01 —0.359 = 0.720 - 
1.182E—01 5.323E—02 2.221 0.026 1.01 

-1.195E—03 2.949E—-04  -4.053 0.000 1.00 

-T.9STE-—O1l  4.584E—01 —1.736 0.083 1.01 
0.354 

—6.996E—01 5.004E—01 —1.398 0.162 - 
1.406E—01 5.346E— 02 2.631 0.009 1.03 
4.887E—02 1.516E-01 0.322 0.747 1.30 

-1.237E-0O1 4.101E—-02 —-3.016 0.003 1.75 
7.956E—02 8.054E—02 0.988 0.323 1.89 

—1.164E+00 4.754E-01 —2449 0.014 1.08 
0.274 
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b) Model 2 
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Fig. 3 Model-based relationships of the number of wolf-killed dogs to the estimated number of wolves occupying the territory and prey abun- 
dances in Finnish wolf territories, 2016-2020. Models are shown in Table | 


fit (AAIC =3.2), largely because white-tail density was 
the primary determinant of ungulate density. The model 
with total ungulate biomass was not weaker than a model 
with a total density of ungulates (model 5 in Table 1, 
AAIC=4.4). 

The model where densities of all three ungulates were 
treated as independent variables (model 6 in Table 1) 
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accounted for 47.7% of the variation in the number of dogs 
killed by wolves. The model with white-tailed deer den- 
sity as the only variable for prey abundance accounted for 
38.4% of the variation, but fit to the results was remarkably 
better (AAIC = 28.0) than that of the model with densities 
of all three ungulate species. 
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In a model where densities of all three ungulate spe- 
cies were entered as independent variables (model 6 in 
Table 1), the number of wolf-killed dogs was related sig- 
nificantly to the density of white-tailed deer but not to the 
density of moose and roe deer. 


Discussion 


The rate of fatal attacks by wolves on domestic dogs was 
related to the abundance of wild ungulate prey in Finland’s 
wolf territories. To what extent the risk of attacks is related 
to ungulate abundance remains, however, somewhat ques- 
tionable because the number of dogs that are accessible 
to wolves within wolf territories is not explicitly known. 
The abundance of dogs is probably correlated with human 
densities that are far higher in southwestern than in east- 
ern Finland. Hunting with dogs is a key issue here because 
almost all casualties were connected to hunting with dogs. 
Unleashed dogs are used in many kinds of hunting in Fin- 
land, most commonly hunting for moose, hare, and grouse. 
For roe deer, dogs are not used that often, and white-tailed 
deer are mostly shot from blinds at agricultural fields and 
feeding sites without using dogs in hunting. 

Although the attacks were highly related to white-tailed 
deer abundance, we had to be cautious with conclusions 
about the risk because the data on the number of dogs that 
can be attacked by wolves are incomplete. The density of 
human residences is higher in wolf territories located in 
regions where white-tailed deer is abundant, and conse- 
quently, dogs are very likely less abundant outside the distri- 
bution range of this deer. Unfortunately, the official register 
for damages launched in 2010 did not provide records of 
whether the attack had taken place in house yards or in hunt- 
ing situations. To have an idea of how often wolves kill dogs 
just in house yards in a situation where they are dependent on 
low-density moose population for food (Gade-Jérgensen and 
Stagegaard 2000; Kojola and Kuittinen 2002), we explored 
data recorded in North Karelia during 1998-2009 by Finn- 
ish Wildlife Agency (J. Kuittinen, unpublished). In 30 cases 
(26%), wolves had killed a dog in a house yard. When related 
to the number of wolf packs in North Karelia (cf. Kojola 
et al. 2014), the mean number of dogs killed in house yards 
was 0.61. The annual mean cases/packs ratio that took into 
account all cases in southwestern Finland during 2017-2022, 
regardless of whether they had occurred in house-yards 
or hunting situations, was 0.39 (I. Kojola unpublished). 
Although practices in the performance of hunting may differ 
regionally, we do not have a reason to suggest that the lower 
use of hunting dogs could explain why so much fewer dogs 
were killed by wolves within southwestern wolf territories 
where white-tailed deer was abundant because the use of 
hunting dogs differ by target species but probably much less 


Page7of9 26 


by the region. Landscapes in southwestern wolf territories 
are most fragmented by farmland and highly trafficked roads 
(Kaartinen et al. 2015). Fragmentation may ease the loca- 
tion of wolves in winter, but on the other hand, the prevail- 
ing dense network of forest roads in eastern Finland which 
wolves commonly used as their movement routes (Gurarie 
et al. 2011) probably promotes the detection of wolf tracks. 

We focused on casualties taken place within wolf territo- 
ries. Vagrant wolves may, however, be more motivated to kill 
dogs than territorial wolves. The proportion of kills outside 
the delineated wolf territories was 42%. We do not have data 
on the proportion of vagrants in autumn, but their proportion 
in wolf populations is 10-20% without a specification for the 
season (Fuller et al. 2003). 

Our results provided evidence that the risk to lose the dog 
to wolves is impacted by an abundant prey, white-tailed deer. 
This pattern could be explained by interference competition 
that increases wolves’ motivation to eliminate dogs within 
their territory. The attacks on dogs take place near terri- 
tory boundaries much more often than expected based on 
the position of GPS-collared territorial wolves (Tikkunen 
and Kojola 2019). Such a spatial pattern fits to an explana- 
tion that intraguild competition is an important motivator 
for wolves to kill dogs. Regardless of the primary reason 
(competition, preying) moose being the only important prey 
could result in a higher risk of attacks. 

A holistic approach in wildlife management that encom- 
passes food chains and species interactions may limit dam- 
ages large carnivores cause to domestic animals (Meriggi 
et al. 1996; Cervasi et al. 2014; Khorozyan et al. 2015; 
Janeiro-Oteroa et al. 2020). Presently high densities of 
ungulates cause substantial damages to forestry and agri- 
culture in many locations in Europe (Valente et al. 2015), 
but at the same time, positive effects of their abundance 
on the viability of large carnivore populations are obvious 
(Chapron et al. 2014). 

Our results provide an example of how a high deer den- 
sity could influence the risk at which wolves attack domestic 
dogs. The negative relationship between the deer density 
and the number of dogs killed by wolves might not hold true 
only in our study system. In years when roe deer and wild 
boar (Sus scrofa) populations are down, the risk of wolves 
attacking dogs in house yards has been shown to be high in 
Estonia (Kojola et al. 2022). 


Management implications 


It is notifiable that high deer density might decrease wolf 
attacks on dogs. The overall applicability might, however, 
remain limited because the key species in our study system, 
white-tailed deer, is an alien species whose present high 
densities are largely owing to winter feeding by hunters. It 
comes expensive to the society at its present density due 
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to damage to forestry and agriculture and being annually a 
partner in > 6000 traffic collisions (Statistics Finland 2022). 

Decreasing the risk of wolf attacks on dogs in regions 
of single-species prey population in eastern Finland pre- 
sumably would be a key in the mitigating conflict between 
wolves and hunters. Higher moose densities which could 
be achieved through decreasing the hunting harvest of 
the moose population would, on the other hand, increase 
the risks of traffic collisions and browsing damages to 
forestry. 
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